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ABSTRACT 

Context. Gravitational lens modeling is presented for the first discovered example of a three-component source for which each com- 
ponent is quadruply imaged. The lens is a massive galaxy member of the cluster CI J0152.7-1357 at z - 0.84 . 

Aims. Taking advantage of this exceptional configuration and of the excellent angular resolution of the HST Advanced Camera for 
Surveys (ACS), we measure the properties of the lens. In particular, the tensing mass estimates of the galaxy are compared to those 
from stellar dynamics and multiwavelength photometry. 

Methods. Several parametric macroscopic models were developed for the lens galaxy, starting from pointlike to extended image mod- 
els. By combining tensing, stellar dynamics, photometry, and spectroscopy, we find an allowed range of values for the redshift of the 
source and the required minimum amount of dark matter enclosed within the disk defined by the Einstein ring of the lens. 
Results. For a lens model in terms of a singular isothermal sphere with external shear, the Einstein radius is found to be 
Re = 9.54 + 0.15 kpc. The external shear points to the cluster's northern mass peak. The unknown redshift of the source is determined 
to be higher than 1.9 and lower than 2.9. Our estimate of the lensing projected total mass inside the Einstein radius, Miq^(R < 9.54 kpc), 
depends on the source distance and lies between 4.6 and 6.2 x 10^^ M©. This result turns out to be compatible with the dynamical 
estimate based on an isothermal model. By considering the constraint on the stellar mass-to-light ratio that comes from the evolution 
of the Fundamental Plane, we can exclude the possibility that at more than 4 cr level the total mass enclosed inside the Einstein ring 
is only luminous matter. Moreover, the photometric- stellar mass measurement within the Einstein radius gives a minimum value of 
50% (1 cr) for the dark-to-total matter fraction. 

Conclusions. The lensing analysis has allowed us to investigate the distribution of mass of the deflector, also providing some in- 
teresting indications on scales that are larger (cluster) and smaller (substructure) than the Einstein radius of the lens galaxy. The 
combination of diff'erent diagnostics has proved to be essential in determining quantities that otherwise would have not been directly 
measurable (with only the currently available data): the redshift of the source and the amount of dark matter in the lens. 

Key words, cosmology: observations - galaxies: high-redshift - galaxies: elliptical and lenticular, cD - dark matter - gravitational 
lensing - galaxies: kinematics and dynamics 



1. Introduction 

Accurate measurements of the mass present in galaxies in the 
forms of dark and visible matter define the empirical framework 
for the theory of galaxy formation and evolution. The knowledge 
of the galaxy mass function is essential for testing galaxy forma- 
tion models (e.g., Hernquist & Springel [20031 De Lucia et al. 
120061) . Among the estimators of total mass in elliptical galaxies, 
stellar dynamics is particularly valuable in the local Universe 
(e.g., Saglia et al. fT992l Gerhard et al. I2QQT1) . while gravita- 
tional lensing oflTers the best direct diagnostics at high redshift 
(e.g., Kochanek 1995; Impey et al. 1998) . Recently, spectral en- 
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ergy distribution (SED) fitting methods have achieved a level 
of precision sufl&ciently high to give reliable estimates of the 
mass contained in the form of stars (e.g., Fontana et al. 120041 
Rocca-Volmerange et al. .2004.; Saracco et al. 120041 Grillo et al. 
2008a). By combining (e.g., Trott & Webster l2QQ2l Treu et al. 
2006; Koopmans et al. 2006) or by comparing (e.g., Drory et al. 
2004; Ferreras et al. 2005 ; Rettura et al. 2006 ) diflferent mass es- 
timators, the internal structure of galaxies can be investigated. In 
particular, the total density distribution and the fraction of mass 
in the form of dark matter can be determined. 

Gravitational lensing has proved to be a unique tool to 
measure the (total) projected mass of galaxies and clusters of 
galaxies on radial scales from kiloparsecs to megaparsecs (e.g., 
Kochanek |1995t Broadhurst et al. .2005 ), and in particular, to 



2 



C. Grillo et aL: A twelve-image gravitational lens system in the z - 0.84 cluster CI J0152.7-1357 




study in detail the relationship between dark and luminous mat- 
ter for these systems. 

In this respect, strong lensing models of multiply imaged 
sources are especially important because the large number of 
observational constraints can lead to robust and detailed mass 
models. The most complex system known so far is a ten-image 
radio gravitational lens (Sykes et al. 1998; Nair 1998 ). 

CI JO 152.7- 1357 is a rich, irregular, and X-ray luminous dis- 
tant i z - .84) cluster (see Delia Ceca et al. [20001 Demarco 
et al. "2005). Its inner regions exhibit spectacular strong grav- 
itational lensing features in the form of multiple images and 
arcs. Strong lensing of background galaxies is also detected on 
smaller scales, around some individual cluster members. A new 
and particularly interesting example is studied in this paper. 

The paper is organized as follows. In Sect. 2, the obser- 
vations of the galaxy cluster CI JO 152.7- 1357 are described. 
In Sect. 3, we address the strong lensing analysis of the sys- 
tem presented in Fig. [T] We report the derived information on 
the redshift of the multiply imaged source in Sect. 4. Then, in 
Sect. 5, different stellar and total mass measurements of the 
lens galaxy are compared. Finally, in Sect. 6 we summarize 
the results obtained in this paper. Throughout this work we 
assume the following values for the cosmological parameters: 
Hq = 70 km s"^ Mpc"\ = 0.3, and Qa = 0.7; in this model 
V corresponds to a linear size of 7.57 kpc at the lens plane. 



2. Observations 

CI J0152.7-1357 was discovered in the ROSAT Dee p Clus ter 
Survey (RDCS; Rosati et al. 1998; Delia Ceca et al. l2QQQl) as 
an extended source with a double core structure. Spectroscopy 
of six galaxies (Ebeling et al. .2000 ) confirmed the cluster and 
gave a redshift of z - 0.83. The X-ray properties of CI JO 15 2.7 - 
1357 have been studied in greater detail by BeppoSAX (Delia 
Ceca et al. 2000), XMM -Newton, and Chandra (Maughan et al. 
120031 Huo et al. 12004b . The X-ray observations were used to 
derive the X-ray luminosity of the cluster and the temperature 
{kT ^ 6 keV, with typical errors of 20%) and metallicity of the 
intra-cluster medium (ICM). The same data provided some ev- 
idence of a possible merger of the two main subclumps, one to 
the northeast (hereafter northern subcluster) and the other to the 
southwest (hereafter southern subcluster). These two structures 
are at a projected distance of 1.6' (corresponding to 730 kpc 
at the cluster redshift.) 

The cluster was observed with the Wide Field Channel of 
the ACS in November and December 2002. It was imaged in the 
F625W (r), F775W (/), and F850LP (z) bandpasses as part of 
a Guaranteed Time Observation program (proposal 9290). The 
observations were done in a 2x2 mosaic pattern allowing for a 
^ 50'' overlap between pointings, with integrated exposure of ^ 
4800 s (see Jee et al. [20n5 ; Blakeslee et al. [2006b . By using AC^ 
observations, a weak lensing analysis provided a detailed mass 
map of the cluster (Jee et al. 2005). The mass reconstruction 
showed complicated substructure in spatial agreement with the 
two peaks of the ICM traced by the X-ray emission. Moreover, 
the X-ray (2.4+^^ x lO^^^©) and weak lensing [(2.1 ± 0.3) x 
10^"^ Mq] estimates of the total mass enclosed within small radii 
(^50'') turned out to be consistent. 

Further multi-band optical and near-IR imaging observa- 
tions of CI JO 152. 7-1357 were performed to select targets for 
spectroscopy (see Demarco et al. 2005; J0rgensen et al. 2005). 
Optical photometry in the B-, V-, and /-bands was obtained 
with the Low Resolution Imaging Spectrometer (LRIS) at the W. 




Fig. 1. Left: Color image of a small portion (~ l.T x ~ 0.7') 
of the northern subcluster created using the riz HST/ACS filters. 
The green circle on the top right, at a projected distance of 488 
kpc and in direction 26.4° NE with respect to the centroid of the 
northern mass clump, which corresponds to two central galax- 
ies (indicated by the red cross on the bottom left), shows the 
position of the strong lensing system considered in this paper. 
Right: Zoom-in view of a 4'' x 4'' field around the lens galaxy at 
zi = 0.82. 



M. Keck Observatory. The LRIS images cover a region of 4.9' x 
6.5'. The cluster was observed in the V-, R-, and / -bands with 
the FOcal Reducer and low dispersion Spectrograph (FORSl) 
at the VLT. The FORSl images cover a region of 6.8'x 6.8'. 
Imaging was also done in the three r'-, /'-, and -bands with the 
Gemini Multi-Object Spectrograph on Gemini North (GMOS- 
N). The GMOS-N images cover a region of 5.5' x 5.5'. The clus- 
ter was imaged in the near-infrared (/- and ^^-bands) with SofI, 
on the NTT at the Cerro La Silla Observatory. The SofI images 
cover a region of 4.9' x 4.9'. 

CI JO 152. 7-1357 was at the center of an extensive spectro- 
scopic campaign carried out with both FORSl and F0RS2 at the 
VLT (Demarco et al. .2005) . A total of 1 1 masks for multi-object 
spectroscopy were used, covering the wavelength range 4000- 
10000 A, with a total exposure time of ^ 29.1 hr. More than 
200 redshifts were measured and 102 galaxies were confirmed as 
cluster members. Spectroscopic observations in the wavelength 
range 5000-10000 A were also performed with one mask of the 
GMOS-N, with a total exposure time of 21.7 hr (J0rgensen et 
al. 120051) . Redshifts and central velocity dispersions were mea- 
sured for 41 galaxies, 29 of which were identified as cluster 
members. From these analyses it was found that the irregular 
distribution of the cluster members follows the extended X-ray 
emission. Furthermore, the mean redshift of the cluster was mea- 
sured z = 0.837 ± 0.001, and the values of velocity dispersion 
of ^ 919 and ^ 737 km s"^ were estimated for the northern and 
southern subclusters, respectively. 

In conclusion. X-ray, weak lensing, and dynamical studies 
all agree in finding in this system that light traces total mass, 
which is concentrated in two main peaks. The inner parts of these 
subclusters coincide with the two regions where strong gravita- 
tional lensing mostly occurs. 

Finally, the temperature and the mass [M(< 65") = (2.1 ± 
0.7) X 10^"^ Mq] of the cluster were also measured from Sunyaev- 
Zel'dovich eff'ect observations, obtained with the interferometers 
of the Berkeley-Illinois-Maryland Association (BIMA ; Joy et al. 
2001 ). A good agreement with the results inferred from the X- 
ray analyses was found. 

In Fig. [T] we show a color-composite image of a small re- 
gion in the northern subcluster obtained by combining the three 
HST/ACS filters. The strong lensing system considered in this 
work is indicated by a green circle, on the left, and shown in de- 
tail, on the right of Fig.[T] The cluster lens galaxy is located at a 
projected distance of 488 kpc from the centroid of the northern 
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Table 1. The lens galaxy. Table 4. The best-fit parameters for the four-image models. 



R.A. Dec. zi e P.A. tab 

(J2000) (J2000) n (mag) {") 

01:52:42.43 -13:56:17.96 0.82 0.08 51.7 22.84 0.28 



Model 












x' 


dof^ 




n 


n 


n 




n 






SIE'^ 


1.91 






0.486 


28.3 


14.3 


3 


SIE 


2.55 


0.15 


-0.17 


0.699 


28.1 


1.77 


1 


SIS+ES b 


1.22 






0.198 


29.1 


6.70 


3 


SIS+ES 


1.23 


0.06 


0.05 


0.185 


28.6 


1.00 


1 



^ Ellipticity or external shear values, depending on the model. 
^ The lens center is fixed to the galaxy center. 
^ Number of degrees of freedom. 



and the position angle 6e for an SIE; the shear y and its posi- 
tion angle 6y for an SIS+ES (see Keeton 2001). Initially, each 
model is analysed by fixing the lens center (xi^ X2;) to the mea- 
sured galaxy position; in a second stage, the lens coordinates are 
also taken as free parameters. Varying these parameters and the 
position of the source, we minimize the chi-square function 

where x^^^^ is the position vector of the i-th observed image (see 
Table|2]), is the corresponding position predicted by the model, 
and cr^ is the position uncertainty. This last quantity is fixed to 
one image pixel (0.05'0 because the centroids of the extended 
images are not well defined. 

The best-fit parameters are shown in Table IH These mod- 
els have diff'erent properties. The SIS is not able to reproduce 
the correct number of images, and thus the parameters of this 
model are not registered. The SIE, with the center coordinates 
as free parameters, has a low;^^, but the high value of the lens 
ellipticity, compared to that measured for the galaxy (see Table 
[B, suggests that this model is inadequate. The least value is 
achieved by an SIS+ES, setting the lens center as a free param- 
eter; the second best model is the SIS+ES with fixed center (the 
higher is partially compensated for by the higher number of 
degrees of freedom, dof). In Fig. [2l we plot the source and im- 
age planes with caustics and critical curves for the last model of 
Table m The positions of the four images, typical of a cusp con- 
figuration, are there compared with the positions predicted by 
the model. We also show the Fermat potential (see Schneider et 
al. 1992 ), its stationary points, corresponding to the image loca- 
tions (two minima and two saddle points with, respectively, the 
same and the opposite parity of the source), and the predicted 
level of distorsion of the images from a round source. 

We notice that the high values found for the external shear 
are plausible in a massive cluster, as is the case of CI JO 152.7- 
1357. For instance, a lens cluster with an Eistein radius of 200 
kpc causes an external shear of 0.2 at a distance of 500 kpc from 
the cluster center (a simplified SIS model is assumed here to 
model the lens). We will see in Sect. 5 that the values of the 
external shear determined in this section are fully consistent with 
the mass estimates of the northern subcluster measured through 
studies of X-ray emission (Huo et al. 2004 ), weak tensing (Jee et 
al. l2QQ5b , and cluster member dynamics (Demarco et al. l2QQ5b . 
Moreover, it is interesting to remark that the position angle of the 
ellipticity for the SIE and that of the shear for the SIS+ES differ 
by about one degree. This common direction points towards the 
closest cluster mass peak, coinciding with the northern clump 
center of the cluster shown in Fig.[Tl 



Table 2. Photometry of the four images. 
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1.07 25.99 


C 
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0.97 25.87 



^ With respect to the galaxy center. 
Table 3. Astrometry of the sub-images. 
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-0.28 


1.04 
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D2 
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XI 
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0.72 
-0.93 


0.66 
-1.10 


0.90 
-0.77 


-0.34 
-0.98 


-0.24 
-0.93 


-0.54 
-0.88 



^ With respect to the galaxy center. 



mass clump (which corresponds to two central galaxies), where 
the most multiple images and arcs produced by CI JO 152.7- 1357 
were detected. This geometrical configuration will justify the 
need for an external shear component in the tensing models of 
the next section. From its spectral and photometrical properties, 
the lens can be safely classified as a bright elliptical galaxy lo- 
cated at a redshift zi of 0.82 and with a stellar central velocity 
dispersion (Xo of 239 + 21 km s"^ (see J0rgensen et al. 2005 ). 
Despite the outstanding observational eff'orts that will be de- 
scribed in Sect. 4.2, the redshift of the quadruply-imaged source 
Zs was not secured. It was only possible to give an upper limit: 
Zs < 2.9. The measured photometric properties of the tensing 
system are summarized in Tables [T]|2l andO more details on the 
properties of this system will be given below. 

3. Gravitational lensing models 

At first, we model the source as a single point-like object, which 
is tensed approximately in the four images of Table [S] Then, 
we employ a three-component source, trying to reproduce the 
twelve images of Table [3] (each of the previous four images, 
hence also the source, is considered here as a triple object). 
Finally, we study models with extended sources, in order to com- 
pare the model-predicted luminosity distribution of the images 
directly with the observations. 

3.1. Four images 

For the lens we consider three diff'erent models: a singular 
isothermal sphere (SIS), a singular isothermal ellipsoid (SIE), 
and a singular isothermal sphere with external shear (SIS+ES). 
An SIS is characterized by three parameters: the Einstein angle 
Oe and the two source coordinates (^1,^2). The remaining mod- 
els involve two additional parameters: the ellipticity e = I -b/a 



C. Grillo et aL: A twelve-image gravitational lens system in the z - 0.84 cluster CI J0152.7-1357 



1.0 r 



0.5 



Source plane 



Image plane 




-1.0 L 

-1.0 -0.5 0.0 0.5 

(arcsec) 



■ m 
A 






B/ 






\ D 




\ n 





Table 5. The best-fit parameters for the twelve-image models. 



-10 12 

(arcsec) 



Image plane 



Image plane 



b 




Fig. 2. Best four-image SIS+ES model. Top left: Source plane 
with caustics. The predicted source position is indicated by the 
cross. The arrows display the directions of the northern cluster 
clump (^Nc), of the external shear predicted by the model {6y), 
of the galaxy ellipticity position angle (P.A.), and of the southern 
cluster clump (^sc)- Top right: Image plane with critical curves. 
The observed (squares) and predicted (crosses) image positions 
are shown. Bottom left: Contour levels of the Fermat potential. 
The images are two minima (A, C) and two saddle points (B, D). 
Bottom right: Predicted deformation of the four images obtained 
by locating a small round source in the position illustrated on the 
top left panel. 



3.2. Twelve Images 

We start here from a triple source (labelled by an index running 
from one to three) lensed into twelve images angularly close to 
the measured positions of Table O Considering the results of the 
previous subsection, a three-component source, for which each 
component is imaged four times, is the most natural assumption 
on the source structure. Three of the observed images, B/, are 
so close to each other that they cannot be deblended; hence, we 
just measure their bright center. The position uncertainty of each 
sub-image is chosen to be one image pixel, as before, except 
for B/ which is taken five times larger, by analyzing the peculiar 
luminosity distribution of the pixels. We model the lens as an 
SIE and an SIS+ES. 

The;^^ minimization results are summarized in Table [51 The 
best;^^, 23.3, is achieved by an SIS+ES, a model consisting only 
of three parameters. If the lens center is also taken as a free pa- 
rameter, this reduces the;^^ by 1 .7, but lowers by two the number 
of dof. The source and image planes with caustics and critical 
curves for the third model of Table [5] are plotted in Fig.[3l 

We note that the value of the Einstein angle of the SIS is 
very similar to that of the best four-image model and that the 
values of magnitude and orientation of the external shear give 
good evidence about the northern subcluster mass peak, as al- 
ready mentioned in the last subsection. 
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^ Ellipticity or external shear values, depending on the model. 
^ The lens center is fixed to the galaxy center. 
^ Number of degrees of freedom. 
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Fig. 3. Best twelve-image SIS+ES model. Left: Source plane 
with caustics. The predicted components of the source are 
shown. Right: Image plane with critical curves. The observed 
(squares) and predicted (triangles, crosses, and diamonds) im- 
ages are shown. 



Then, we determine the statistical uncertainties on the 
Einstein angle and on the external shear, starting from two hun- 
dred minimizations on simulated data sets. In particular, we 
extract image positions from Gaussian distributions centered on 
the measured values and with standard deviations equal to the 
position uncertainties reported previously, and for each set we 
search for the best-fit parameters. 

The parameter distributions are shown in Eig.|4]and the re- 
sults are presented in Table[6l The marginal density functions are 
approximately Gaussian and the joint density functions reveal 
correlations (for definitions see Cowan 1998): the correlation 
coefficients are ro^^y = -0.64, re^^e^ = -0.03, and ry^e^ = 0.19. 
The origin of the significant anti-correlation of the value of the 
Einstein angle and the magnitude of the external shear will be ex- 
plained in Sect. 5.1. We notice here that, at zi = 0.82, an Einstein 
angle of (1.26 + 0.02) corresponds to an Einstein radius (Re) 
of (9.54 + 0.15) kpc. Moreover, we remark that the low uncer- 
tainty values on the parameters suggest that the measurements 
are accurate and robust. 

The model just presented is valuable in describing the lens 
properties with only three parameters, although it cannot repro- 
duce the complex image configuration very accurately. In order 
to improve the agreement between the observed and the recon- 
structed image geometry we have tried more refined models, 
without gaining any real improvement. For instance, some el- 
lipticity is inserted in the SIS+ES model, but, as expected from 
the known degeneracy between external shear and ellipticity (see 
Witt et al. |1997K we cannot find a lower value of the;^^. 

We interpret the external shear as due to the mass distribu- 
tion on the cluster scale. So, instead of employing an external 
shear component, we decide to model both the northern subclus- 
ter (the most relevant mass clump for this study) and the galaxy 
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Table 7. The best-fit parameters for the extended-image model. 



0.15 0.16 0.17 0.18 0.19 0.20 0.21 



Fig. 4. Results of the;^^ minimizations of 200 Monte-Carlo gen- 
erated data sets. Top: Marginal distributions of the model param- 
eters with their approximate normal distributions. Bottom: The 
joint density distributions show the parameter correlations. 

Table 6. The best-fit parameters with the relative marginalized 
errors for the twelve-image SIS+ES model. 
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1.26 + 0.02 0.179 + 0.011 26.4 + 1.1 



as SISs. The description of the subcluster in terms of an SIS is 
only a first-order approximation, but it is adequate to represent 
the mass distribution at distances from the center large compared 
to the Einstein radius of the clump (below we will prove that this 
statement is valid in our system). The two Einstein angles are the 
only needed parameters (hence, the number of dof is here 16). 
In fact, the first center is set equal to the average position of the 
two brightest galaxies of the northern subcluster (see Fig.[T]), the 
second one to the center of the lens galaxy. We measure a 
of 25.5, and Einstein angles for the subcluster (^Ec) ^^d for the 
galaxy (^Eg) of 19.6'' (7?Ec = 148 kpc) and 1.07'' (/?Eg = 8.10 
kpc), respectively. This latter value will be used in the following 
as a conservative lower limit of the galaxy Einstein angle. 

Finally, as far as the source properties are concerned, the 
predicted three-component source has a linear extent of about 
0.09", equivalent to approximately 0.7 kpc at a hypothetical red- 
shifts of 2.5. In addition, the total absolute magnification of the 
sub-images {ji = I/zaI + lyUel + \i^c\ + ly^ol, where jii is the mag- 
nification of the /-th image, with a positive or negative sign if 
the image is, respectively, a minimum or a saddle point of the 
Fermat potential) is 27.9, a reasonable value in the vicinity of 
critical lines. 

Taking into consideration the irregular and compact shape 
and the blue color of the source, it is very likely that this object 
is a high-redshift star-forming region. 

3.3. Extended Images 

As a final step, we employ an extended parametric approach. 
For this purpose, we have developed an ad hoc algorithm, which 
reconstructs the properties of both a lens and a source by com- 
paring the observed and the predicted luminous intensity of an 
array of pixels on the image plane. Before presenting the results, 
we describe this technique shortly, since it is not as well-known 
as the point-like parametric (Keeton 2001) or the non-parametric 
(Saha & Williams 1 1997i; Koopmans , .20051) approaches. 
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^ The lens center is fixed to the galaxy center. 



In principle, the method is straightforward, because it ex- 
ploits only the ray-tracing equation, yix) = x - a(x ), and the 
surface brightness conservation in lensing, I(x ) = P{y(x )) (for 
more details see Schneider et al. 119921) . First, the image and 
the source planes are divided into sub-pixels. In particular, each 
CCD pixel of the image plane is splitted into 16 sub-pixels. 
Furthermore, in order to deal with large magnified images, the 
pixel size on the source plane is fixed to 1/256 the pixel size of 
the CCD. Then, starting values for the parameters of the mod- 
els that describe the lens (/?/) and the luminosity distribution 
of the source (ps ) are assigned. According to the tentative lens 
model, each point of the image grid is ray-traced into the source 
plane and is associated with the closest point of the source grid. 
Finally, the images of the lensed source are obtained by giv- 
ing to the image sub-pixels the values of the source sub-pixels 
matched in the previous step. The eff'ect of the point spread func- 
tion (PSF), I(x) = (PSF >^ P){y(x)), and the effect of the sub- 
pixels (anti-aliasing) are taken into account. The best parameters 
for the adopted lens and source models are found by minimizing 
the following chi- square function 



(/obs(^)-/(^)) 



0-] 



A^pix 



(2) 



where A/pix is the total number of pixels of the observation con- 
sidered in the modeling, /obs(-'^) and I{x) are, respectively, the 
intensity observed and predicted by the model in the CCD pixel 
located at jc, and cri is the standard deviation of the intensity eval- 
uated on a blank field of the CCD (i.e., the noise). The method 
has been tested on artificial but plausible lens systems and has 
provided good results for the values of the reconstructed param- 
eters. 

Following the best models discussed above, we represent the 
lens as an SIS+ES and, on first approximation, we model the 
source as the sum of three gaussian functions with independent 
parameters. 

The best-fit parameters and model resulting from the applica- 
tion of our algorithm are shown in Table [7] and Fig. [51 The prop- 
erties of the lens are essentially the same as those obtained pre- 
viously; on the other hand, we note slight variations in the values 
of the source parameters. In fact, the values of the Einstein an- 
gle and of the external shear are consistent with those of the best 
point-like model (see Table [6]). The largest separation among the 
centers of the three gaussians now gives a predicted linear extent 
for the source of about 0.18", which corresponds to approxi- 
mately 1.4 kpc at a redshift of 2.5. The total absolute magnifi- 
cation, i.e. the integrated intensity of the images divided by that 
of the source, is on the order of 14.9. Table [8] summarizes the 
signed (positive or negative if the nature of the two Fermat sta- 
tionary points of each pair is, respectively, the same or not) flux 
ratios of the images A (/a//c), B (/b//c), and D (/d//c), with 
respect to the most magnified image C, and the cusp ratio (7?cusp) 
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Fig. 5. Best extended-image SIS+ES model. Top left: The ob- 
served HST/ACS r field (5'' across). Top right: The same image 
after the subtraction of an elliptical model fitted on the luminos- 
ity profile of the lens galaxy. Some residuals of the galaxy lumi- 
nosity are still visible in the inner region. Bottom left: The best 
reconstruction of the lensed system. Bottom right: The residuals 
after the subtraction of the model predicted (third panel) from 
the observed (second panel) images. The central residuals due to 
the lens galaxy subtraction have been masked here. 

Table 8. Signed flux and cusp ratios. 





/a//c 


/b//c 


/d//c 


(^cusp)bCD 


Obs. 


0.79 + 0.11 


-0.54 + 0.08 


-0.60 + 0.08 


0.07 + 0.06 


Mod. 


0.41 


-0.59 


-0.56 


0.07 



for the closest triplet of images B, C, and D. This last quantity, 
which is introduced by Keeton et al. l2QQ3l can be written as 



(^cusp)bCD 



l + f 

1 Jc 


-h 


/d| 
/c| 


1 + 


|/b| 

|/c| 


-h 


/d| 
/c| 



(3) 



In Table [8] we report the observed values and their uncertainties, 
assuming a 10% error on each flux, and the model-predicted val- 
ues of the previously defined quantities. For the three images B, 
C, and D, we measure an opening angle ^ of 93.6° and a dimen- 
sionless separation d/O^ of 1.20 (for definitions see Keeton et al. 
[20031) . 

The image positions taken into consideration in the pointlike 
models are supposed to coincide approximately with the centers 
of the images of the extended models. Nonetheless, the former 
analyses consider essentially only the regions where the images 
are maximally amplified by lensing, whereas the latter studies 
probe also the smaller amplification regime. Therefore, we do 
not have any obvious explanation for the diff'erent predicted sizes 
of the source; while, the diff'erence between the values of the 
total magnification found here and in the previous paragraph can 
be easily understood. 



For an ideal cusp catastrophe, the cusp ratio should exactly 
vanish (Schneider et al. 11992b , but, since real lenses are not 
ideal cusps, this general property is expected to hold only ap- 
proximately, depending on the opening angle and the dimension- 
less separation of the triplet considered (for a detailed study see 
Keeton et al. l2003 D. Surprisingly, most of the observed lens sys- 
tems show a significant violation of the cusp relation (known 
also as flux ratio "anomaly"; e.g., see Kent & Falco 1988; Sluse 
et al. 12003b . This fact can be explained by assuming the presence 
of small-scale structure in the lens galaxy, on a scale comparable 
to the separation between the images (Mao & Schneider 1 19981) . 
Although we measure a cusp ratio which is consistent with zero 
(and thus, given the measured values of the opening angle and 
of the dimensionless separation, the cusp relation is here satis- 
fied), a possible clue on the existence of small perturbations in 
the lens potential might be the significant inconsistency between 
the observed and the predicted flux ratio for image A (see Table 
[8] and the last panel of Fig. [5]). 

4. The source redshift 

A precise value for the redshift of the lensed source is not avail- 
able yet. Its measurement is very difficult because of the faint- 
ness of the object and its proximity to the bright lens galaxy (see 
Table [21). However, some information about this redshift can be 
extracted from different diagnostics, as explained below. 

4.1. Gravitational lensing and stellar dynamics combined 
The Einstein angle for an SIS is explicitly given by 

\2 



Dis 
Do 



(4) 



where crsis is the lens "velocity dispersion", c is the speed of 
light, Dis and Dos are the lens-source and the observer- source 
angular diameter distances, respectively. It follows from Eq. Q 
that 6e may be interpreted as a function of both the lens mass 
M (M oc (T^j^r^) and the source redshift Zs (Dis, Dos)- A lower 
limit for this redshift can be obtained by combining gravitational 
lensing and stellar dynamics. 

In fact, several tests (Kochanek [T993l [T9941 Treu et al. [20061 
Grillo et al. 2008b) have proved that in elliptical galaxies the 
stellar central velocity dispersion (ctq) is a good estimate of the 
velocity dispersion of a one-component isothermal model ((Xsis). 
A spectroscopic measurement of ctq for the lens galaxy is pre- 
sented by J0rgensen et al. l2Q05l ((Tn = 239 ±21 km s"^). Starting 
from this result and from the value of the galaxy Einstein angle 
of the previous section, in Fig. [51 the velocity dispersions from 
gravitational lensing and from stellar dynamics are found to be 
compatible at a confidence level of 95% if Zs > 2.1 and at 99% 
if > 1.9. 



4.2. Spectroscopic information 

The spectra of objects A, B, and C were obtained with 
VLT/F0RS2, as part of an extensive observational campaign car- 
ried out on the galaxy cluster (see Demarco et al. 12005 1 for 
more details on the spectroscopic analysis of other objects in CI 
JO 152. 7-1357). An unambiguous estimate of the source redshift 
was not possible from the spectra of images B and C, due to the 
contamination from the overwhelming lens galaxy light in 0.8'' 
seeing conditions. The spectrum of image A was taken in three 
different exposures, with a total integration time of 5230 s, with 
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Fig. 6. The source redshift obtained from gravitational lensing 
and stellar dynamics. The Einstein angle, from lensing analy- 
sis, and the velocity dispersion of the galaxy, from spectroscopic 
measurements, combine to give a lower limit for the redshift of 
the source. The dashed and dotted lines represent 1 and 2 cr error 
bars, respectively. The vertical long-dashed line shows the spec- 
troscopic upper limit of the source redshift discussed in Sect. 4.2 





Fig. 7. Spectroscopic measurements of object A. Top: Position of 
the r'slit. Bottom: Composite 2D sky- subtracted spectrum ob- 
tained with VLT/F0RS2 in 5230 s. The arrow shows the Lyman 
alpha drop at 4730 A. 



the 300V grism, using the same mask geometry and the same slit 
characteristics. The position of the slit and the combined 2D sky- 
subtracted spectrum are shown in Fig. [71 Although the signal is 
low, no significant emission feature is observed, the continuum 
is visible down to 4730 A, however the cross-correlation with 
template spectra does not yield satisfactory results. If we inter- 
pret this limit as the Lyman alpha drop, blueward to the Lyman 
alpha emission line, then we can infer that the redshift of the 
source (zs) is lower than 2.9. 

5. Mass measurements 

In this section we first determine the lens galaxy total (lumi- 
nous+dark) mass, based on the previously discussed lensing 
models and on the published (J0rgensen et al. 12005b stellar 
velocity dispersion measurement. Next, we measure the lens 
galaxy stellar mass, based on multiwavelength photometry spec- 
tral template fitting. This allows us to infer a lower limit for the 
dark matter fraction enclosed inside the Einstein radius of the 
deflector. 



Table 9. Luminous and dark mass measurements of the lens 
galaxy, without the cluster contribution. The ranges represent the 
intervals at 68% confidence level. 





Mien 


Mdyn 


Mphot 


/dm 




(10^1 Mo) 


(10^1 Mo) 


(10^1 Mo) 




/?< S.lOkpc 


3.9-5.2 


2.8 - 4.0 


0.6 - 0.9 


0.44 - 0.99 


R < 9.54 kpc 


4.6 - 6.2 


3.3-4.7 


0.7-1.1 


0.50 - 0.99 



Galaxy mass estimate inside 




SIS+SIS mass estimates 



IQn _ \ m,JR< 9.54 kpc) 



1.0 



MjJR< 9.54 kpc) 



1.5 2.0 2.5 3.0 3.5 

Source redshift (zj 



McusteXR ^ 488 kpc) 



Mcakny^ClusreXR^ 9.54 kpc) 1 

Mc,u,JR< 9.54 kpc) 



4.0 



1.0 



Mc,„,„,(^< 9.54 kpc) 

1.5 2.0 2.5 3.0 3.5 

Source redshift (z) 



4.0 



Fig. 8. The galaxy and cluster masses. Left: The lensing 
[Mien(/? < 9.54 kpc)] and dynamical [Mdyn(^ < 9.54 kpc)] mea- 
surements of the galaxy projected total mass inside a 9.54 kpc 
radius as a function of the source redshift. The dashed and dot- 
ted lines represent 1 and 2 cr error bars, respectively. Right: 
The galaxy (dashed), the cluster (dotted), and the total (dashed- 
dotted) projected mass inside a 9.54 kpc radius centered on the 
galaxy. The long-dashed line shows the projected mass of the 
cluster inside a circle with a radius of 488 kpc, that is the dis- 
tance between the cluster and the galaxy centers. 



5.1. Lensing estimates 

A measurement of the value of the Einstein radius of a lensing 
system can be directly translated into a mass estimate for the de- 
flector. In particular, the Einstein radius is defined as that radius 
inside which the mean surface density of a lens is equal to the 
critical surface density (Zcr) of the system (for more details see 
Schneider et al. 119921) . So, it follows that the projected lensing 
mass (Mien) enclosed within Re equals 



Mien(< ^e) = ^crnRl 



(5) 



In our specific case, Scr can assume a range of values depending 
on the source redshift. This is the main source of uncertainty in 
our lensing mass estimates, more relevant than that associated 
to detailed modeling of the cluster. From the results of Sect. 3.2 
and 4, we have that the projected total mass of the lens galaxy 
without the cluster contribution Mien(^ < 9.54 kpc) is well mea- 
sured and its value is included between 4.6 and 6.2 xlO^^ Mq, 
at a 68% CL. This information is shown in Table [9] and on the 
left of Fig. [U 

The lower value of Req, found in Sect. 3.2 when considering 
the cluster mass contribution, suggests that the galaxy must be 
less massive than expected from the SIS-hES model. It is inter- 
esting to measure the projected total mass inside Re = 9.54 kpc, 
predicted by the SIS-hES and SIS-hSIS models. In the former 
case, we use exactly Eq. ([5]); in the latter case, the mass is the 
sum of the cluster and of the galaxy masses inside the same cir- 
cle (see Fig. [SI on the right). The diff'erence between the values of 
the total mass from the two models is less than 0.3%. Therefore, 
the lensing mass estimate is found to be very reliable and in- 
dependent of the model details. To sum up, the galaxy and the 
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cluster masses inside the galaxy Einstein radius may have differ- 
ent weights, but their sum is nearly constant. This explains the 
anti-correlation noted above between the Einstein radius and the 
external shear, which is interpreted as the cluster contribution. 
Moreover, the projected mass of the cluster inside the circle of 
radius 7? = 488 kpc, the projected distance between the centers 
of the northern subcluster and the lens galaxy, agrees well with 
the value of about 3 x 10^^ Mq from analyses in X-rays (Huo 
et al. l2QQ4l) , weak lensing (Jee et al.l2QQ5]) , and dynamics of the 
cluster galaxies (Demarco et al. 120051) . 

5.2. Dynamical estimates 



Table 10. Multi-band photometry of the lens galaxy. 



Several studies (e.g., Rusin et al. 120031 Koopmans et al. 12006b 
have shown that the total (luminous+dark) density distribution 
of elliptical galaxies is homologous and well described by a 
(isothermal) profile. A similar conclusion was also reached by 
stellar dynamical studies (out to Rq) of nearby galaxies, based on 
the application of self-consistent equilibrium dynamical models 
that incorporate the picture of galaxy formation by collisionless 
collapse (Bertin & StiavelH 1993 ). 

We consider here mass measurements projected along the 
line of sight. The projected mass (M) enclosed inside a certain 
radius (R) for an SIS is 



M(< R) = 



(6) 



where osis has been introduced in Eq. (|4l), and G is the universal 
gravitational constant. As mentioned in Sect. 4.1, a good estima- 
tor of (Tsis is found to be ctq, and this latter quantity has been 
measured for the galaxy studied in this paper. Starting from Eq. 
© and the value of ctq, we evaluate the dynamical projected total 
mass, Mdyn(< R), within the two radii considered in the lensing 
analysis of the previous section. The results are presented in the 
third column of Tableland on the left of Fig. [8l 

The lensing and dynamical mass estimates are consistent and 
the lensing measurements do not turn out to be more accurate 
than the dynamical ones only because a rough estimate of the 
source redshift is available. 



5.3. Photometric estimates 

An estimate of the photometric- stellar mass (Mphot) of a galaxy 
can be derived by comparing the observed SED of the galaxy 
with a set of composite stellar population (CSP) templates, com- 
puted with stellar population models. In addition to the stel- 
lar mass, this method also allows the age and the star forma- 
tion history (SFH) of the galaxy to be investigated. It is well- 
known that the derived stellar mass depends on the adopted ini- 
tial mass function (IMF) and only weakly on the assumed model 
of dust extinction and metallicity evolution. Moreover, it has 
been shown (Rettura et al. I2006D that the photometric- stellar 
mass does not exhibit statistically significant discrepancies when 
evaluated with diff'erent stellar population models (e.g., Bruzual 
& Chariot 2003 vs. Maraston 2005). Finally, it must be noted 
that accurate stellar mass measurements require unbiased galaxy 
SEDs, which translates into accurate PSF-matched photometry. 

We derive the photometric- stellar mass of the lens galaxy 
through this multi-wavelength matched aperture photometry 
method. We use Bruzual & Chariot's (2003) templates at solar 
metallicity, assuming a Salpeter (1955 ) time-independent IMF 
and a delayed exponential SFH. First, we smooth all images (r, 
/, z from HST/ACS and /, Ks from NTT/Sofl) to the worst PSF 







e< 1.07" 


e< 126" 


HST/ACS 




23.688 + 0.014 


23.441 + 0.013 


HST/ACS 


lAB 


22.415 + 0.005 


22.171 +0.005 


HST/ACS 


Zab 


21.755 + 0.004 


21.511+0.004 


NTT/SofI 


Jab 


21.219 + 0.010 


20.990 + 0.009 


NTT/SofI 


Ksab 


20.426 + 0.011 


20.183 + 0.011 
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Observed Wovelength A (/xm) 
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Fig. 9. SED of the lens galaxy at z/ = 0.82. The circles with the 
error bars represent, from left to right, the observed flux densities 
measured inside the 1.07''aperture in the HST/ACS (r, /, z) and 
NTT/Sofl (/, Ks) passbands. The illustrated best-fit has been built 
with Bruzual & Chariot 2003 models. On the bottom, the best-fit 
values of the age (t), the characteristic time of the SFH (r), and 
the mass (Mphot) of the galaxy are given. 



and then we measure the galaxy magnitudes inside the apertures 
suggested by the lensing analysis (see Sect. 3.2). The results are 
summarized in Table [TOl In Fig. [9] we plot the observed SED of 
the galaxy and the best-fit Bruzual & Chariot 2003 model for 
the first aperture. The best-fit intervals, at the 68% confidence 
level, of the stellar mass inside the two apertures are shown in 
the fourth column of Tabled 

Starting from the lensing and dynamical measurements of 
the total mass and from the photometric measurement of the stel- 
lar mass, we have determined the minimum amount of dark mat- 
ter (/dm). The radii relative to which /dm is computed are about 
four times larger than the eff'ective radius {Rq) of the galaxy (see 
Table [B. In order to obtain the dark matter fraction, we assume 
the "maximum light" hypothesis: inside each aperture the stellar 
mass is divided by the lower total mass estimate (in this particu- 
lar case, the latter value comes from the dynamical study). This 
gives an upper limit to the stellar mass fraction, from which a 
lower limit to the dark mass fraction can be inferred. In the last 
column of Table [9] the measured values of /dm are presented. 

The total luminosity of the system inside the 9.54 kpc 
Einstein radius is (3.1 ± 0.3) x lO^^Lo^g, implying an average 
mass-to-light ratio M/Lg = (13.0 ± 2.6) Mq L~^^ (considering 
the lower dynamical estimate of the total mass). Under appro- 
priate assumptions (e.g., see Treu et al. 120011 [200 5) the evolu- 
tion of the intercept of the Fundamental Plane with redshift can 
be related to the evolution of the average eff'ective mass-to-light 
ratio. The SLACS and LSD Surveys (Treu et al. 2006) have es- 
tablished that, in the redshift range from to 1 , the eff'ective stel- 
lar mass-to-light ratio evolves as dlog(M / Lb) / dz = -0.76, with 
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an rms scatter of 0.11. If the evolution of the effective mass-to- 
light ratio is equal to the evolution of the stellar mass-to-light 
ratio, we can use this result to infer M^/Lb of our lens galaxy, 
assuming that log(MJLB)z = log(MJLB)o + Alog(M/Lfi). The 
first term on the right-hand side of the previous equation can 
be measured for local E/SO galaxies; e.g., using the data from 
Gerhard et al. |2001j (M,/Lfi)o = (7.9 ± 2.3) MoL'^ Hence, 
we find that MJLb (z = 0.82) = (1.9 ± 0.7) MqL-]^. We note 
that the value M^/Lb = 13.0 MqL~^^ required to explain the 
projected mass enclosed inside the Einstein radius solely by lu- 
minous matter is inconsistent with the above independently de- 
rived value at more than 4 cr level. This is further evidence of the 
presence of a significant dark matter component inside the lens 
Einstein circle. Moreover, from the estimate of the photometric- 
stellar mass we obtain for M^/Lb a value between 2.3 and 3.5, 
which is consistent with the value derived from the evolution of 
the Fundamental Plane. 

We mention that our findings on the fraction of the dark mat- 
ter component are compatible with the results from other diff'er- 
ent studies on the amount of dark matter in elliptical galaxies 
(e.g., Saglia et al. 1992; Treu & Koopmans"2004'; Ferreras et al. 
120051 Treu et al. 2006 ). We also note that our relatively high val- 
ues of /dm support the picture that the most massive ellipticals 
be dark matter dominated in their outer regions. 



6. Conclusions 

In this paper we have presented a strong lensing analysis of 
an elliptical galaxy, member of the cluster CI JO 152. 7-1357 
(z - 0.84). By means of HST/ACS deep observations, a three- 
component source highly magnified into twelve images was dis- 
covered. For this system we have discussed several paramet- 
ric macroscopic models: first point- like, then extended image 
models have been developed. The Einstein radius of a singular 
isothermal sphere was found to be Re = 9.54 ± 0.15 kpc, and 
the value of the projected mass inside this radius was shown to 
be a function of the unknown source redshift Zs- By combining 
lensing with galaxy dynamics and by studying the source spec- 
trum, we have obtained for this redshift lower and upper bounds 
of 1.9 and 2.9, respectively. These limits have allowed us to 
get accurate and robust total mass measurements of the galaxy: 
Mien(/? < 9.54 kpc) = (4.6 - 6.2) x 10^^ Mq. Furthermore an ex- 
ternal shear component was proved to be indicative of the north- 
ern mass distribution of the cluster and the predicted value of the 
cluster mass has been shown to be in agreement with the X-ray, 
weak lensing, and cluster member dynamical analyses. Then, we 
have measured the lens total and luminous mass from stellar dy- 
namics, and from optical and near-IR photometry. From these 
results we have estimated a lower limit of 50% (1 cr) for the 
dark matter fraction enclosed inside the Einstein radius of the 
galaxy. The presence of a significant dark matter component is 
also confirmed by comparing the value of the mass-to-light ra- 
tio measured in our lens with the stellar one predicted from the 
evolution of the Fundamental Plane. 

Using non-parametric methods (see Bradac et al. il2005l) , the 
complex image configuration of the lensing system might be bet- 
ter reproduced, but possibly with only a small gain of insight into 
the lens properties. 
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